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Hilbert space are analysed. Assuming that the corresponding classical dynamics takes place on 
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define the Monge distance between two arbitrary density matrices. The Monge metric has a simple 
semiclassical interpretation and induces a non-trivial geometry. Among all pure states the distance 
from the maximally mixed state p*, proportional to the identity matrix, admits the largest value 
for the coherent states, while the delocalized 'chaotic' states are close to p*. This contrasts the 
geometry induced by the standard (trace, Hilbert-Schmidt or Bures) metrics, for which the distance 
from p* is the same for all pure states. We discuss possible physical consequences including unitary 
time evolution and the process of decoherence. We introduce also a simplified Monge metric, defined 
in the space of pure quantum states and more suitable for numerical computation. 
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I. INTRODUCTION 



Consider two quantum states described by the density matrices pi and p2- What is their distance in the space 
of quantum states? One should not expect a unique, canonical answer for this question. On the contrary, several 
possible distances can be defined, related to different metrics in this space. As usual, each solution possesses some 
advantages and some drawbacks; each might be useful for different purposes. 

Perhaps the simplest possible answer is given by the norm of the difference. The trace norm leads to the trace 
distance 



Ar(pi,P2) =tr -P2)2 (1.1) 

used by Hillery to describe the non-classical properties of quantum states and by Englert [|| to measure the 
distinguishability of mixed states. In a similar way the Frobenius norm results in the Hilbert-Schmidt distance 



Dhs{puP2) = v/tr[(pi ~ P2)2] (1.2) 

often used in quantum optics 

Another approach based on the idea of puriGcation of a mixed quantum state leads to the Burcs distance . An 
explicit formula for the Bures distance was found by Hiibner |9| 



,(pi,P2) = \/2(l~tr[(p}/W}/^)i/2]), (1.3) 



and various properties of this distance are a subject of a considerable interest (see [|10|-|16[). It was shown by Braunstein 
and Caves [[17| that for neighbouring density matrices the Bures distance is proportional to the statistical distance 



introduced by Wootters |18| in the context of measurements which optimally resolve neighbouring quantum states. 

Note that for pure states pi = \ipi){ipi\ and p2 = \^2){^2\ we can easily calculate the above standard distances, 
namely 



A.(l¥'i),l¥'2)) =2v/T^, (1.4) 

Dhs{\^i),\V2))^vW^, (1-5) 

and 

DBures{\^l),\^2)) ^ ^ 2 {1 - ^) , (1.6) 

where the transition probability p = \{(pi\ip2)\'^ = cos^(S/2). The angle S equals to the Fubini-Study distance 
Dpsi^i, ^2) in the space of pure states, and for iV = 2 it is just the angle between the corresponding points of the 
Bloch sphere |jl^. The Fubini-Study metric defined by 

Dpsifi, V2) = arccos(2p — 1) = 2 arccos-^, (1.7) 

corresponds to the geodesic distance in the complex projective space (see e.g. [^) and for infinitesimally small values 
of p becomes proportional to any of standard distances. 

In a recent paper [ pO[ we introduced the Monge metric Dm in the space of density operators belonging to an 
infinite-dimensional separable Hilbert space Ti.. The Monge metric fulfils the following semiclassical property: the 
distance between two harmonic oscillator (Glauber) coherent states |ai) and \a2) localised at points ai and 02 of the 
classical phase space = C is equal to the Euclidean distance d between these points 

DM(|cki), |a2)) = d(ai, a2) . (1.8) 
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In the semiclassical regime this condition is rather natural, since the quasi-probability distribution of a quantum 
state tends to be strongly localised in the vicinity of the corresponding classical point. A motivation to study such 
a distance stems from the search for quantum Lyapunov exponent, where a link between distances in the Hilbert 
space and in the classical phase space is required |21 ] . Our construction was based on the Husimi representation of a 
quantum state p given by 

Fp(«):=i(a|p|a), (1.9) 

TT 

for a e C. The Monge distance Dm between two arbitrary quantum states was defined as the Monge-Kantorovich 
distance between the corresponding Husimi distributions [ pO| . 

Although the Monge-Kantorovich distance is not easy to calculate for two or more dimensional problems, 



it satisfies the semiclassical property (1.8), crucial in our approach. On the other hand, one could not use 
for this purposes any 'simpler' distances between the Husimi distributions, like e.g. Li or L2 metrics, be- 
cause the semiclassical property does not hold in these cases. Moreover, this property is not fulfilled for 
any of the standard distances in the space of density matrices (trace, Hilbert-Schmidt or Bures distances). 
Consider two arbitrary pure quantum states |(^i) and [(^2) G ^ and the corresponding density operators 
Pi = and p2 = \V2){V2\- If the states are orthogonal, the standard distances between them do not de- 

pend on their localisation in the phase space. For example the Hilbert-Schmidt and the Bures distances between two 
different Fock states \n) and |m) are equal to v^, and the trace distance is equal to 2. Although the state |1) is localised 
in the phase space much closer to the state |2) then to |100), this fact is not reflected by any of the standard distances. 
Clear ly, t he same concerns a nonlinear function of the Hilbert-Schmidt distance, which satisfy the semiclassical condi- 
tion ( |l.8| ) and was recently introduced in |2^]. On the other hand, the Monge distance is capable to reveal the phase 
^ace structure of the quantum states, since Dm(|to), \n)) = \arn ~ an\, where Uk — y/TT(^^^{2k + 1)72^*^+^ ~ \/k (see 

In this paper we propose an analogous construction for a classical compact phase space and the corresponding finite- 
dimensional Hilbert spaces Hn. In particular we discuss the N — {2j + 1)— dimensional Hilbert spaces generated by 
the angular momentum operator J. In the classical limit the quantum number j tends to infinity and the classical 
dynamics takes place on the sphere 5*^. However, the name Monge metric on the sphere should not be interpreted 
verbally: the metric is defined in the space of density matrices, while the connection with the sphere is obtained via 
the SU{2) vector coherent states, used in the construction to represent a quantum state by its generalised Husimi 
distribution. In general, the Monge distance in the space of quantum states can be defined with respect to an arbitrary 
classical phase space fi. 

This paper is organised as follows. In Sect. II we review some properties of pure and mixed quantum states in 
a finite-dimensional Hilbert space. In Sect. Ill we recall the definition of the Monge metric based on the Glauber 
coherent states and extend this construction to an arbitrary set of (generalised) coherent states. We analyse basic 
properties of such defined metric and its relation to other distances in the space of density operators. The case where 
the classical phase space is isomorphic with the sphere 5^, corresponding to the 5/7(2) coherent states, is considered 
in Sect. IV. We compute the Monge distance between certain pure and mixed states, and compare the results with 
other distances (trace, Hilbert-Schmidt, and Bures). In particular, we give the formulae for the Monge distance 
between two coherent states (for arbitrary j) and between two arbitrary mixed states for j — 1/2. In the latter case 
the geometry induced by the Monge distance coincides with the standard geometry of the Bloch sphere induced by 
the Hilbert-Schmidt (or the trace) distance. However, in the higher dimensions both geometries differ considerably. 
Potential physical consequences of our approach are discussed in Sect. V. In Sect. VI we introduce a simplified version 
of the Monge metric, defined only in the space of pure quantum states, but better suited for numerical computation. 
Finally, some concluding remarks are provided in Sect. VII. 



II. SPACE OF MIXED QUANTUM STATES 



A. Topological properties 

Let us consider a pure quantum state \tp) belonging to an A^— dimensional Hilbert space Hn- It may be described by 
a normalised vector in T^at, or by the density matrix = Such a state fulfills the purity condition: p'^ = p^. 

The manifold V, containing all pure states, is homeomorphic with the complex projective space CP^~^. This space is 
2(A^ — 1)— dimensional. In the simplest case N = 2, the two-dimensional space CP^ corresponds to the Bloch sphere. 
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To generalise the notion of pure states one introduces the concept of mixed quantum states. They are represented 
hy N X N positive Hermitian matrices p, which satisfy the trace condition trp = 1. Any density matrix may be 
diagonaUzed and represented by 



p:^VEV'', (2.1) 

where V is unitary, while a diagonal matrix of eigenvalues E contains only non-negative entries: Ei > 0; i — 1, . . . , N . 
For each pure state all entries of E are equal to zero, but one equal to unity. Due to the trace condition X^i^i — ^■ 
It means that the set of all such matrices E forms an {N — 1)— dimensional simplex Sn in R^. Let i? be a diagonal 
unitary matrix. Since 

p = VEV^ = VBEBW^ , (2.2) 

therefore the matrix V is determined up to N arbitrary phases entering B. On the other hand, the matrix E is 
defined up to a permutation of its entries. The form of the set of such permutations depends on the character of the 
degeneracy of the spectrum of p. 



a) 

N = 2 



Kii 



(10) 



Hi) 



(01) 



(10) 




(100) (010) (100) (iio) 



c) 

N = 4 (0001) 




FIG. 1. (A'^ — 1)— dimensional simplex iSjv of diagonal density matrices of size A'^ and its antisymmetric part Sn for a) A'^ = 2, 
b) TV = 3, and c) TV = 4. The simplex Sn, enlarged at the right hand side, can be decomposed into 2^ parts. TV numbers in 
brackets denote coordinates in the original TV— dimensional space of eigenvalues. Corners of iSjv represent pure states (density 
matrices of rank one), edges - matrices of rank two, faces - matrices of rank two. Bold lines (grey faces) symbolise boundary 
of 5jv. 
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Representation ( ^.2[ ) makes the description of some topological properties of the (A''^ — 1)— dimensional space A4 
easier [^4[^. We introduce the following notation. We write = {Si)'^ = [U{1)]^ for the TV— dimensional torus. 
Identifying points of Sn which have the same coordinates (but ordered in a different way) we get an asymmetric 
simplex Sn- Equivalently, one can divide Sn into A'^! identical simplexes and take any of them. The asymmetric 
simplex Sn can be decomposed in the following natural way 

Sn= U ^fei,...,fc„, (2.3) 

fciH hk„=N 

where n — 1, . . . , N denotes the number of different coordinates of a given point of Sn, fci the number of occurrences 
of the largest coordinate, ^2 the second largest, etc. Observe that Kki,...,k„ is homeomorphic with the set Gn, where 
Gi is a single point, G2 a half-closed interval, G3 an open triangle with one edge but without corners, and generally, 
Gn is an (n — 1)— dimensional simplex with one (n — 2)— dimensional hyperface without boundary (the latter is 
homeomorphic with an (n — 2)— dimensional open simplex). There are N ordered eigenvalues: Ei > E2 > ■ ■ ■ > En, 
and N —1 independent relation operators 'larger or equal', which makes all together 2^~^ different possibilities. Thus, 
Sn consists of 2^"^ parts, out of which (^^Z^) parts are homeomorphic with Gm, when m ranges from 1 to N. The 

decomposition of the asymmetric simplex Sn is illustrated in Fig. 1 for the simplest cases = 2, 3, and 4. 

Let us denote the part of the space M related to the spectrum in Kk-^^...^k„ {n- different eigenvalues, the largest 
eigenvalue has fci multiplicity, the second largest k2 etc.) by Mki,...,kn- ^ mixed state p with this kind of the spectrum 
remains invariant under arbitrary unitary rotations performed in each of the fc.;— dimensional subspaces of degeneracy. 
Therefore the unitary matrix B has a block diagonal structure with n blocks of size equal to fci , . . . , fc„ and 

Mk,^...,k„ ^ [U{N)l{U{ki) X • • • X C/(fc„))] X Gn , (2.4) 
where fci + • • • + fc„ = and fc^ > for i = 1, . . . , n. Thus Ai has the structure 

y Mk,^...,k„ ~ U [U{N)l{U{ki) X • • • X C/(fc„))] X Gn , (2.5) 

fciH Vkn=N fciH \-kn=N 

where the sum ranges over all partitions of A^. The group of rotation matrices B equivalent to F = U{ki) x U{k2) x 
• • • X U{kn) is called the stability group of U{N). 

For A^ = 2 we have Mi^i ^ [U{2)/T^] x G2 ~ S*^ x G2 and M2 ~ {p*}, so the space M. has the topology of a 
two-dimensional ball - the Bloch sphere and its interior. This case and also cases A^ = 3,4 are analysed in detail in 
Tab. 1. 



A^ 


Label 


Decomposition 


Subspace 


Part of the 
asymmetric simplex 


Topological Structure 


Dimension 

D = D1+D2 


1 


Ml 


1 


El 


point 


[C/(l)/C/(l)]xGi = {p.} 


= + 


2 


Mn 


1 + 1 


El > E2 


line with left edge 


[U{2)/T'] X G2 


3 = 2 + 1 


M2 


2 


El = E2 


right edge 


[C/(2)/C/(2)J X Gi ^ {p4 


= + 


3 


Mm 


1 + 1 + 1 


El > E2 > Es 


triangle with base 
without corners 


[C/(3)/r3] X Ga 


8 = 6 + 2 


M12 


1 + 2 


El > E2 = Es 


edges with 
lower corners 


[C/(3)/(C/(2) X T)\ X G2 


5 = 4+1 


M21 


2 + 1 


El = E2 > Es 


Ms 


3 


El = E2 = Es 


upper corner 


[C/(3)/C/(3)J X Gi - {p4 


= + 


4 


Mini 


1+1+1+1 


El > E2 > Es > E4 


interior of tetrahedron 
with bottom face 


[C/(4)/r4] X G4 


15 = 12 + 3 


M112 


1 + 1 + 2 


El > E2 > Es — E4 


faces without side edges 


[C/(4)/({/(2) X T^)] X G3 


12 = 10 + 2 


M121 


1 + 2 + 1 


El > E2 — Es > E4 


M211 


2 + 1 + 1 


El = E2 > Es > E4 


Mrs 


1 + 3 


El > E2 ~ Es — E4 


edges with lower corners 


[C/(4)/(C/(3) X T}\ X G2 


7 = 6 + 1 


Msi 


3 + 1 


El = E2 = Es > E4 


M22 


2 + 2 


El = E2 > Es = E4 


[C/(4)/(t/(2) X U{2}}\ X G2 


9 = 8 + 1 


Mi 


4 


El = E2 = Es = E4 


upper corner 


[C/(4)/[/(4)J X Gi = {p4 


= + 



Table 1. Topological structure of the space of mixed quantum states for a fixed number of levels A^. The group 
of unitary matrices of size A^ is denoted by U{N), the unit circle (one-dimensional torus ~ U{1)) by T, while G„ 
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stands for a part of an (n — 1)— dimensional asymmetric simplex defined in the text. Dimension D of the component 
-^fei,...,fc„ equals Di + D^, where Di denotes the dimension of the quotient space U{N)/T, while D2 = n — 1 is the 
dimension of the part of the eigenvalues simplex homeomorphic with G„ . 

Note that the part Mi^...^i represents generic, non-degenerate spectrum. In this case all elements of the spectrum 
of p are different and the stability group H is equivalent to an iV-torus 

A^i,...,i ^ [C/W/T^] X Gat. (2.6) 

Above representation of generic states enables us to define a product measure in the space of mixed quantum 
states. For this end, one can take the uniform (Haar) measure on U {N) and a certain measure on the simplex Sjq 
1^,0. The coordinates of a point on the simplex may be generated by squared moduli of components of a 
random orthogonal (unitary) matrix [ p9[ . 

The other 2^~^ — 1 parts of M. represent various kinds of degeneracy and have measure zero. The number of 
non-homeomorphic parts is equal to the number P{N) of different representations of the number TV as the sum of 
positive natural numbers. Thus P{N) gives the number of different topological structures present in the space M. 
For iV = 1, 2, . . . , 10 the number P{N) is equal to 1, 2, 3, 5, 7, 11, 15, 22, 30 and 42, while for larger N there is described 

by the asymptotic Hardy-Ramanujan formula |3^, P{N) ^ 

In the extreme case of iV-fold degeneracy, Ei = 1/A^, the subspace Mn ^ [U{N)/{U{N) x T°)] x Gi ^ Gi, so 
it degenerates to a single point. This distinguishes the maximally mixed state := I/A^, which will play a crucial 
role in subsequent considerations. For the manifold of pure states n = 2 and ki = l,k2 = N — 1 (since Ei = 1, 
E2 = ■ ■ ■ = En ^ 0) and so V ^ [U{N)/{U{N - 1) x f/(l))] x (1, 0, . . . , 0) ~ CP^-\ In the case iV = 2 it can be 
identified with the Bloch sphere S'^. 

On the other hand, it is well known that Ai itself has a structure of a simplex with the boundary contained in the 
hypersurface detp = 0, with rank 1 matrices (pure states -V) as 'corners', rank 2 as 'edges', etc., and with the point 
'in the middle' (see [^,^ for a formal statement and |Q for a nice intuitive discussion). 



Let us mention in passing that the quotient space appearing in (2.4) 

U{N) 



C/(fci) X [/(fca) X • • • X [/(fc„) 

is called a Hag manifold, and in a special case 



(2.7) 



Gr(k,N) := , , ^^^^ ^ (2.8) 

^ ' ' U{k) X U{N-k) ^ ' 

a Grassman manifold. For a fuller discussion of the topological structure of Ai (especially for = 4) we refer the 
reader to |2l . 



B. Metric properties 

The density matrix of a pure state = |'0)('0| may be represented in a suitable basis by a matrix with the first 
element equal to unity and all others equal to zero. Due to this simple form it is straightforward to compute the 
standard distances between p^ and directly from the definitions recalled in Sect. 1. Results do depend on the 
dimension TV, but are independent on the pure state namely 



Ar(p^,P*) = 2- , Dhs{P4,,P*) = - , DBures{p-4,,P*) ^ \h- ■ (2.9) 



In the sense of the trace, the Hilbert-Schmidt, or the Bures metric the 2(iV — 1)— dimensional space of pure states V 
may be therefore considered as a part of the {N'^ — 1)— dimensional sphere centred at p* of the radius r depending 
on N and on the metric used. From the point of view of these standard metrics, no pure state on V is distinguished; 
all of them are equivalent. It is easy to show that the distance of any mixed state from p^ is smaller than r, in the 
sense of each of the standard metrics. Thus the space of mixed states A4 lays inside the sphere ~^ embedded in 
although, as discussed above, its topology (for iV > 2) is much more complicated than the topology of the 
{N'^ — 1)— dimensional disk. 
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The degree of mixture of any state may be measured, e.g., by the von Neumann entropy S = — Trplnp = 
— '^f^iEiln{Ei). It varies from zero (pure states) to In(A^) (the maximally mixed state p*). Let us briefly dis- 
cuss a simple kicked dynamics, generated by a Hamiltonian represented by a Hermitian matrix H of size A'^. It maps 
a state p into 

p' = e'"pe-'" , (2.10) 

where the kicking period is set to unity. 

Such a unitary quantum map does not chang e the eigenvalues of p, so the von Neumann entropy is conserved. In 



particular, any pure state is mapped by (2.10) into a pure state. Any mixed state p, which commutes with H, is 
not affected by this dynamics. Assume the Hamiltionian H to be generic, in the sense that its N eigenvalues are 
different. Then its invariant states form an [N — 1)— dimensional subspace 2h C M, topologically equivalent to Sn- 
In the generic case of non-degenerate Hamiltonian it contains only N pure states: the eigenstates of H. Note that 
the invariant subspace Xh always contains p* . 

Moreover, the standard distances between two states are conserved under the action of a unitary dynamics, i.e. 



D,ipi,p2)^D,ip[,p'^), (2.11) 



where denotes one of the distances: Dtr, Djjs or Dsures- Therefore, the unitary dynamics given by ( ^.10 ) can 



be considered as a generalised rotation in the {N'^ — 1)— dimensional space Ai, around the (N — 1)— dimensional 
'hyperaxis' Th, which is topologically equivalent to the simplex S^~^. In the simplest case, iV = 2, it is just a 
standard rotation of the Bloch ball around the axis determined by H. For example, if = aJz, where Jz is the third 
component of the angular momentum operator J, it is just the rotation by angle a around the z axis joining both 
poles of the Bloch sphere. The set Xh of states invariant with respect to this dynamics consists of all states diagonal 
in the basis of J^: the mixed states with diag(p) = {a, 1 — a} (a g (0, 1)) and two pure states, |l/2, 1/2) for a = 1, 
and |l/2,-l/2) for a = 0. 



III. MONGE DISTANCE BETWEEN QUANTUM STATES 

A. Monge transport problem and the Monge-Kantorovich distance 

The original Monge problem, formulated in 1781 |34j, emerged from studying the most efficient way of transporting 
soil [||: 

Split two equally large volumes into inEnitely small particles and then associate them with each other so that the 
sum of products of these paths of the particles over the volume is least. Along which paths must the particles be 
transported and what is the smallest transportation cost? 

Consider two probability densities Qi and Q2 defined in an open set fl C i?", i.e., Qi > and Jq Qi{x)d"x = 1 for 
i = 1, 2. Let Vi and V2, determined by Qi, describe the initial and the final location of 'soil': Vi = {{x,y) G ft x i?+ : 
< y < Qi{x)}. The integral Jy d^xdy is equal to the unity due to normalisation of Qi. Consider one-to-one 
maps T : — > 57 which generate volume preserving transformations Vi into V2, i.e.. 



Qi{x)^Q2{Tx)\T'{x)\ (3.1) 

for all a; € f2, where T'{x) denotes the Jacobian of the map T at point x. We shall look for a transformation giving 
the minimal displacement integral and define the Monge distance p5| , ^ 

^Af(Qi,Q2) :=inf / \x-T{x)\Q^{x)d^x, (3.2) 

where the infimum is taken over all T as above. If the optimal transformation Tm exists, it is called a Monge plan. 
Note that in this formulation of the problem the 'vertical' component of the soil movement is neglected. The problem 
of existence of such a transformation was solved by Sudakov [0, who proved that a Monge plan exists for Qi,Q2 
smooth enough (see also p^). The above definition can be extended to an arbitrary metric space (fi, d) endowed with 
a Borel measure to. In this case one should put d{x, T{x)) instead of \x — T{x)\ and dm{x) instead of d^x in formula 
(3.2), and take the infimum over all one-to-one and continuous T : — > O fulfilling J^Qidm = Jrp-it^) Q2dm for 
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each Borel set A C fl. In fact we can also measure the Monge distance between arbitrary two probabiHty measures 
in a metric space (fi, d). For fijV - probability measures on (f2, d) we put 

Da/Im, J^) := inf / d{x,T{x))dii{x) , (3.3) 
Jn 

where the infimum is taken over all one-to-one and continuous T : ^ fl such that fi{A) = ^{T^^IA)) for each 
Borel set A C il. To avoid the problem of the existence of a Monge plan Kantorovich introduced in 40s 

the 'weak' version of the original Monge's mass allocation problem and proved his famous variational principle (see 
Proposition 1). For this and other interesting generalisations of the Monge problem consult the monographs by Rachev 
and Riischendorf |3^j4l[] . 

In some cases one can find the Monge distance analytically. For the one-dimensional case, 51 = M, the Monge 
distance can be expressed explicitly with the help of distribution functions Fi{x) = J_^ Qi{t)dt, i = 1,2. Salvcmini 
obtained the following solution of the problem p3|| 

/ + 00 
\F^{x) - F2{x)\dx . (3.4) 
-oo 

S ever al two-dimensional problems with some kind of symmetry can be reduced to one-dimensional problems, solved 



by (3^). In the general case one can estimate the Monge distance numerically relying on algorithms of solving 



the transport problem, often discussed in handbooks of linear programming |4£ 





FIG. 2. Estimation of the Monge distance between two overlapping distributions Q\ and Q2 (a) by the distance between the 
reduced distributions Qi{x) — Qiix) — V(x) (b), where V{x) denotes the overlap. 



According to definition ( |3.2| ) taking an arbitrary map T which fulfils (^j) we obtain an upper bound for the Monge 
distance Dm- Another two methods of estimating the Monge distance are valid for a compact metric space {ft,d) 
equipped with a finite measure m. The first method may be used to obtain lower bounds for Dm- It is based on the 
proposition proved by Kantorovich in 40s (see also |^,^,^). 

Proposition 1. {variational formula for the Monge-Kantorovich metric) 



Dm{Qi,Q2) = max 



f{x){Qi{x) - Q2{x))dm{x) 



(3.5) 
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where the supremum is taken over all / fulfilling the condition \ f{x) — f{y)\ < d(x, y) for all x,y £ (weak contrac- 
tions). 

To obtain another upper bound for Dm we may apply the following simple estimate: 
Proposition 2. 

DMiQuQ2) < iA/2)DLAQuQ2), (3.6) 

where Dli{Qi,Q2) = \Qi{x) — Q2{x)\dm{x) and A — diam(ri). 

The intuitive explanation of this fact is the following. Let v be the volume of the 'overlap' of the probability 
distributions Qi and (52, i-e., = V{x)dm{x), where V — minjQi, Q2}- Then Dm{Qi, Q2) < (1 ~ v)/S.^ because 
the number (1 — v) represents the part of the distribution to be moved and the largest possible classical distance on VL 
is smaller than or equal to A. Moreover, Dj^^ (Qi, Q2) — 2(1 — w), which proves the assertion. Although Fig. 2 presents 
the corresponding picture for the simplest, one-dimensional case, Proposition 2 is valid for an arbitrary metric space. 
For the formal proof see Appendix A. 

B. The Monge distance - harmonic oscillator coherent states 

In we defined a 'classical' distance between two quantum states pi and p2 via the Monge distance between the 
corresponding Husimi distributions Hp^ and Hp^ : 

Dm{pi,P2)-^Dm{Hp,,Hp,) , (3.7) 



where Hp^ are given by formula ( L£ ) . Observe that the family of harmonic oscillator coherent states | a) , parameterised 
by a complex number a, is implicitly present in this definition. 

The Monge distance satisfies the semiclassical property: the distance between any two Glauber coherent states, 
represented by Gaussian Husimi distributions localised at points ai and 02, is equal to the classical distance \ai — 02] 
in the complex plane EQ] . 



C. The Monge distance - general case and basic properties 



The above construction, originally performed for the complex plane with the help of the harmonic oscillator coherent 



states, may be extended to arbitrary generalised coherent states of Perelomov |46 defined on a compact classical phase 
space. Let G be a compact Lie group, G 9 g ^ Rg € Ti. its irreducible unitary representation in the Hilbert space 
H, and T the subgroup of G, which consists of all elements y E G leaving the reference state \k) g V invariant 
(i.e. Ry\K) ~ |k)). Define il = G/T and j?/) = i?,,]/*) for [rj] £ G/T. Note that |1) — \k), where 1 is the group 
G unit. Consider a family of the generalised coherent states 9 ^ jry) e 7^. It satisfies the identity resolution 
In I rf™(f?) = 1/ dimTi, where m is the natural (translation invariant) measure on the Riemannian manifold (17, d) 
normalised by the condition m{Q) — 1, and d is the Riemannian metric on f2. Let us denote by C the manifold of all 
quantum coherent states, isomorphic to fi, and embedded in the space of all pure states V. Note that (77I77) = 1. 

For the SU{k) coherent states the space J7 C is isomorphic to CP'^~^ and m is the natural Riemannian measure 
on Q. Obviously, the dimension of the Hilbert space Hn carrying the representation of the group equals N > k, and 
ii N — k all pure states are SU{k) coherent, and C — V. For example, in the case of SU{2) vector coherent states the 
corresponding classical phase space is the sphere ~ CP^ |46|. In the simplest case N — 2 (or j ~ 1/2) pure states 
are located at the Bloch sphere and are coherent. 

Any quantum state p € M may be represented by a generalised Husimi distribution Hp : fl — > defined by 

Hp{7j) N ■ {7j\p\7^) , (3.8) 

for T] G n, which satisfies 

/ Hp{r,)dm{r,) = 1 . (3.9) 
In particular, for a pure state p = \'d){'d\ {\'d) E V) and 77 e we have 

^^|^)(^|(^) :=^I(^I'7>|'- (3.10) 
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In the sequel we shall assume that for coherent states ji?) G C (i? G 17) the densities tend weakly to the Dirac- 

delta measure in the semiclassical limit, i.e., when the dimension of the Hilbert space carrying the representation 
tends to infinity. 

The Monge distance for the Hilbert space Hn, the classical phase space {ft,d) and the corresponding family of 



generalised coherent states jr/) is then defined by solving the Monge problem in fl, in the full analogy with (3.7) 



Dm{pi,P2) Dm (Hp,, Hp,) . (3.11) 

The distance d{r],T{r])) between the initial point rj and its image T{r]) with respect to the Monge plan has to be 
computed along the geodesic lines on the Riemanni an m anifold ft. 



For compact spaces the semiclassical condition (1.8) for the distance between two coherent states becomes weaker: 



Property A. {semiclassical condition) Let ?7i,772 G ^- Then 

DMi\m)Av2))<divum), (3.i2) 

and 

Dm{\i1i)j\V2)) d{rii,ri2) (in the semiclassical limit), (3.13) 
where d represents the Riemannian distance between two points in f2. 



To demonstrate (3.12) it suffices to take for the transformation T in ( p^ ) the group translation ?72 * (e.g. the 
respective rotation of the sphere in the case of SU{2) coherent states). However, this transformation needs not 
to give the optimal Monge plan. As we shall show in the following section, this is so for the sphere and the SU{2) 
coherent states. On the other hand, in the semiclassical limit (for SU{2) coherent states: j oo), the inequality 



in ( 3.12| ) converts into the equality, in the full analogy to the property ( |l.8| ), valid for the complex plane and the 



harmonic oscillator coherent states. This follows from the fact that the Monge- Kantorovich metric generates the 
weak topology in the space of all probability measures on Q, and the densities tend weakly to the Dirac-delta 

6rfi in the semiclassical limit, for i = 1, 2. 

The Monge distance defined above is invariant under the action of group translations, namely: 

Property B. (invariance) Let a,(3 ^ G and pi,p2 G -M- Then 

Dm{Rp''piRp, Rp''p2Rp) = Dm{pi,P2) ■ (3.14) 



Particularly, 



where 1 is the group G unit. 



7?M(|a),|/3)) =7?M(|r'a),|l)), (3.15) 



The above formulae follow from the definition of the Monge distance (3.3), and from the fact that both the measure 
m and the metric d are translation invariant. 

D. Relation to other distances 



Let pi, P2 G M- We start from recalling the variational formula for the trace distance (see for instance |47 
Proposition 3. [variational formula for Dtr) 

Ar(pi,P2)= sup \iv A{pi - p2)\ , (3.16) 

Il-4||<1 

where the supremum is taken over all Hermitian matrices A such that < 1, and the supremum norm reads 

IIAII = sup{Pa;|l -.xeH^AM < 1} ■ (3.17) 
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Applying Proposition 1 we can prove an analogous formula for the Monge distance 
Proposition 4. {variational formula for Dm) 

-Dm(pi,P2)= max | tr^(pi - p2)| , (3.18) 

L{A)<1 

where the maximum is taken over all Hermitian matrices A with L{A) < 1, and 



L{A) = inf{c : there exists a c-Lipschitzian function f : Q ^ R. such that A= f {v)\v) {v\din{ri)} . (3.19) 

For the proof see Appendix B. This proposition has a simple physical interpretation. It says that the Monge distance 
between two quantum states is equal to the maximal value of the difference between the expectation values (in these 
states) of observables (Hermitian operators) some of whose P-representations are weak contractions. Recently, Rieffel 
p8| considered the class of metrics on state spaces which are generated by Lipschitz seminorms. If il is compact, then 
one can show that the Monge metric Dm belongs to this class. 

From Propositions 3 and 4 we can also easily deduce Proposition 2. Using Proposition 2 and the Holder inequality 
for the trace (see E^) one can prove the following inequalities 



-Dm <Dl,<N- Dhs <N-Dtr, (3.20) 

where A is the diameter of il and N ~ dimHj\f. On the other hand from the fact that the Monge-Kantorovich 
metric generates the weak topology in the space of probability measures on fi, it follows that Dm{pi, P2) implies 
Dhs{pii P2) for every pi,p2 € -M. Thus the Monge metric Dm and the Hilbert-Schmidt metric Dhs generate 
the same topology in the space of mixed states M. 

Let us emphasise here a crucial difference between our 'classical' Monge distance and the standard distances in the 
space of quantum states. Given any two quantum states represented by the density matrices pi and p2, one may 
directly compute the trace, the Hilbert-Schmidt or the Bures distance between them. On the other hand, the classical 
distance is defined by specifying the set of generalised coherent states in the Hilbert space. In other words, one needs 
to choose a classical phase space with respect to which the Monge distance is defined. Take for example two density 
operators of size = 3. The distance Dm{Pi,P2) computed with respect to the SU{2) coherent states and, say, with 
respect to the SU{i) coherent states can be different. The simplest case of the SU{2) coherent states corresponding 
to classical dynamics on the sphere is discussed in the following section. 



IV. MONGE METRIC ON THE SPHERE 

A. Spin coherent states representation 

Let us consider a classical area preserving map on the sphere Q : ^ S'^ and a corresponding quantum map U 
acting in an A— dimensional Hilbert space Ti.N- A link between classical and quantum mechanics can be established 
via a family of spin coherent states {"d^ S H localised at points {■&, (p) of the sphere CP^ — S'^. The vector coherent 
states were introduced by Radcliffe B3] and Arecchi et al. pffl and their various properties are often analysed in 
the literature (see e.g. mJ52]). They are related to the S'[/(2jalgebra of the components of the angular momentum 



operator J = {Jx, Jy, Jz}, and provide an example of the general group theoretic construction of Perelomov |46| (see 
Sect. HI C). 

Let us choose a reference state usually taken as the maximal eigenstate \j,j) of the component acting on 
Hn, N — 2j + 1, j = 1/2, 1, 3/2, .... This state, pointing toward the 'north pole' of the sphere, enjoys the minimal 
uncertainty equal to j. Then, the vector coherent state is defined by the Wigner rotation matrix R^^ip 

1^9,^) = R^Jk) = + (4.1) 

where R^,^, = exp [ii? (cost^Ja: — simpjy)], ~ — iJ^j and 7 = tan(i9/2)e*''', for {-d, ip) £ S*^ (we use the spherical 
coordinates). 

We obtain the coherent states identity resolution in the form 



11 



/ \d,v){d,v\dti{^,^)^I/{2j + l), (4.2) 

where the Riemannian measure dfi{'d, ip) — {sun!} / ATr)d'dd(p does not depend on the quantum number j. 

Expansion of a coherent state in the common eigenbasis of and J^: |j, m), m = — j, . . . , +j (in Hn) reads 



\^,ip)^ ^ sinJ-'"(z9/2)cosJ+™(i5/2)exp(i(j-m)v3) 



2j 

J — TO 



1/2 



(4.3) 



The infinite 'basis' formed in the Hilbert space by the coherent states is overcomplete. Two different SU{2) coherent 
states overlap unle ss th ey are directed into two antipodal points on the sphere. Expanding the coherent states in the 
basis of Hat as in (4.3) we calculate their overlap 



(4.4) 



where S is the angle between t wo v ectors on related to the coherent states and \i}',(p'), and for j = 1/2 it 

equals to the geodesic distance (1.7). Thus, we have 



H\^,^)is^^\{^',ip') = (2j + l)cos4^(S/2). (4.5) 

This formula guarantees that the respective Husimi distribution of an arbitrary spin coherent state tends to the 
Dirac ^-function in the semiclassical limit j —^ oo. 

To calculate the Monge distance between two arbitrary density matrices pi and p2 of size N one uses the N = 
(2j + 1)— dimensional representation of the spin coherent states ip) (to simplify the notation we did not label them 
by the quantum number j). Next, one computes the generalised Husimi representations for both states 



Hp^{^,^):^N-{d,M^,v) (4.6) 

and solves the Monge problem on the sphere for these distributions. Increasing the parameter j (quantum number) 
one may analyse the semiclassical properties of the Monge distance. 

It is sometime useful to use the stereographical projection of the sphere onto the complex plane. The Husimi 
representation of any state p becomes then the function of a complex parameter z. It is easy to see that for any 
pure state lip) £ V the corresponding Husimi representation is given by a polynomial of [N — 1) order: W^{z) = 
^N-i _|_ ^^^p^cjz' = with arbitrary complex coefficients q. This fact provides an alternative explanation of the 
equality V — CP^~^. Thus, every pure state can be uniquely determined by the position of the {N — 1) zeros of 
on the complex plane (or, equivalently, by zeros of on the sphere). Such stellar representation of pure states is 
due to Majorana [ |53t and it found several applications in the investigation of quantum dynamics ||5^-p6[|. In general, 
the zeros of Husimi representation may be degenerated. This is just the case for the coherent states: the Husimi 
function of the state \d, (p) is equal to zero only at the antipodal point and the {N — 1)— fold degeneracy occurs. The 
stellar representation is used in section VI to define a simplified Monge metric in the space of pure quantum states. 



B. Monge distance between some symmetrical states 



Consider two quantum states, whose Husimi distributions are invariant with respect to the horizontal rotation. 
Using Proposition 4 one may found the Monge distance between both states with the help of the Salvemini formula 

Proposition 5. Let pi,p2 G M fulfil Hp.{d,ip) = Hp^{d) for {d,ip) £ S*^. Consider the normalised one- 
dimensional functions hp. : [0,7r] M+ given by hp.{{)) := iJp^ (t?)i sin that satisfy hp.{'&)d'& — 1 for i = 1,2. 
Then 
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Jo 



(4.7) 



where the cumulative distributions read Fp.{'d) = hp.{'4')dip. 

The proof is given in Appendix C. We use this proposition computing the Monge distance between two arbitrary 
eigenstates of the operator Jz (Sect. IV Fl) and the Monge distance of some of these eigenstates from the maximally 
mixed state (Sect. IV F2). 

The maximally mixed state p, = 1/N is represented by the uniform Husimi distribution on the sphere. Thus, 
HpSd,ip) = Hp,{'0) = 1, = ssini?, andFp.(i?) = i(l-cos^'). 

Furthermore, all the eigenstates of , forming the basis | j, to) , possess this symmetry, and according to ( |4.3| ) we 
get 



h\j,^)U,mm = (2J + 1) (^^. \) sin2(^--™)+i(^/2) cos2(^+'")+i(^/2) . 



(4.8) 



The formula (4/7) enables us to compute the Monge distance between them. 
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FIG. 3. Quantum states, for which we calculate Monge distances, are denoted at the Bloch sphere corresponding to j = 1/2. 
Dots at small circles represent the position of zeros of the Husimi function of the corresponding pure states. 



We introduce the notation p+ = |j, j)(j, j|, P- — |j, and pa = ap+ + (1 ^ for a G [0, 1] (for N — 2 

these states are represented in Fig. 3). It follows from Proposition 1 that Hp^ = aHp^ + (1 ^ o)Hp_ , and, consequently, 
DM{p+,Pa) = (1 - a)DM{p+,P-), and DM[Pa,P~) = aDM[p+,P-)- 

In some cases we can reduce the two-dimensional problem to the Salvemini formula, even if it docs not possess 
rotational symmetry. 

Proposition 6. Let pi , p2 G . Define : [0, tt] x [0, 2tt] ^[0,1] by Fi{t,ip) = \ /J Hp^ ip) sin d d-d for t G [0, tt] , 
(f G [0, 27r], and i — 1,2. Assume that 

1. Fi(7r, if) = F2(7r, ip) for aU cp G [0, 27r], 

2. Fi{t, (f) > F2{t, ip) for aU t G [0, tt] and ip G [0, 27r]. 
Then 



DM{pi,P2) = 7r / {Fi{t,p)-F2{t,^))dtdp. 
271- Jo Jo 



(4.9) 



For the proof see Appendix D. We use the above proposition computing both, the Monge distance between two 
arbitrary density matrices for j = 1/2 (Sect. IVC), and the Monge distance between two coherent states for arbitrary 
j (Sect. IV F3). 
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C. Monge distances for j — 1/2 



Let us start from the calculation of the Monge distance between the 'north pole' and a mixed state pa 
parameterised by a G (0,1) (note that ~ pi/2 in this case). Computing the distribution functions we get 
Fp^{-d) = (sin^?9)(2a - l)/4 + (1 - cosi?)/2, while Fp^{-d) = Fp^^^{-d). Elementary integration gives the re- 

sult Dm{p+, Pa) — ~ a)7r/4. Substituting a = 1/2 for p^ or a = for p_ we get two important special cases: 



Dm{p+,P-) = 7r/4; Dm{p+,P*) = Dm{p--,P*) = tt/ 



(4.10) 



These three states p+, and p_ lay on a metric line. This follows from the property of the distribution functions 
visible in Fig. 4. They do not intersect, and therefore the area between and F- is equal to the sum of two 
figures: one enclosed between F+ and F*, and the other one enclosed between and F_. Note, however, that the 
distance Z)m(/3+, P-) ~ 0.785 is much smaller than the classical distance between two poles on the sphere equal to tt. 
Instead of rotating the distribution Hp^ by the angle tt, the optimal Monge plan consists in moving south the 'sand' 
occupying the north pole, along each meridian. The difference between both transformations is so large only in this 
deep quantum regime, for which the distributions are very broad and strongly overlap. As demonstrated below, this 
effect vanishes in the semiclassical regime j — > oo, where the semiclassical property (f]q) is recovered. 



j=l/2 



a) 





FIG. 4. a) Cross-section of the Husimi distribution Hp^i)) for p+ (peaked at the top), for p, (uniform), and for p_ (peaked 
at the bottom) plotted for j = 1/2. b) Monge distance between these states may be represented as the area between graphs of 
the corresponding distribution functions Fp^(-d), Fp^(-d), and Fp_{-d). 

In the case j — 1/2 all pure states are coherent, so the Monge distance from p* is the same for every pure state (as 
illustrated in Fig. 3). Thus, the manifold of pure states (the Bloch sphere) forms, in the sense of the Monge metric, 
the sphere of radius i?i = 7r/8 centred at p*. All mixed states are less localised than coherent and their distance to 
p* is smaller than To see this note that every mixed state can be represented as a vector v in the unit ball. Using 
Proposition 6 and some geometrical considerations one can find the Monge distance between any two mixed states pi 
and p2- Representing them by Pauli matrices a and vectors Vi in the Bloch ball of radius 1/2, namely, pi = p^ + a ■ Vi, 
we obtain [57| 

TT TT 

Dm{pi,P2) = -d{vi,V2) = -\vi -V2\, (4-11) 
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where d denotes the EucHdean metric in R^. Consequently, for j = 1/2 the Monge distance induces the same geometry 
as that of the Bloch ball, as illustrated in Fig. 3. 



D. Monge distances for j = 1 

In an analogous way we treat the case iV = 3. Obtained data 

Ri:^Dm{p+,P*)^Dm{p-,P*) = 3tt/16; Dm{p+, P-) = 3tt/8; (4.12) 

R2 ■■= Dm{p*,P\o)) ^ i/G; Dm{p+,P\o}) = DM{p-,piQ'j) = 3tt/16. (4.13) 

are based on the results derived in Appendix E (see also the next subsection) and visualised in Fig. 5b. Note that 
both triples and {p+,p|o),P-} lay on two different metric lines. Thus, in contrast to the case j — 1/2, 
the two states and p_ are connected by several different metric lines. 




FIG. 5. Schematic map showing the eigenstates \j, m) and the mixed state p, for a) j — 1/2, b) j = 1, c) j = 3/2, d) j = 2, 
and e) j ^ 00. The symbol * labelling dots represents the maximally mixed state p*, while numbers m denote pure states 
I J, m) . Solid lines denote the metric lines, and the accompanying numbers represent the approximate Monge distance between 
the states. 

Now, consider a mixed state pm represented in the canonical basis by a diagonal density matrix pm — diag(a, b, c), 
where a + b + c — 1. Since {p+,p|o),P-} lay on a metric line and their distributions are invariant with respect 
to the horizontal rotation, it is not difficult to calculate the Monge distance Dm{p+, Pm) using Proposition 5. 
The corresponding distribution functions do not cross, and so Dm{p+, Pm) = i>D m{p+, P\o)) + cD m{p+, p ) = 
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11(2 - 6 - 2a). For comparison DHs{p+,Pni) = {I - a)^ + b'^ + {1 - a - fe)^, Dtr{p+,Pm) = \/2(l - a), and 
Dsures — ^72(1 — \/a). This simple example shows that for j — 1 the Monge metrics induces a non-trivial geometry, 
considerably different from geometries generated by any standard metric. 

The Monge distance between any pure state {ip) and the mixed state depends only on the angle x between 
two zeros of the Husimi function located on the sphere. If the zeros are degenerated, x = 0; the state is coherent and 
Rjf, = Ri. The coherent states are as much localised in the phase space, as allowed by the Heisenberg uncertainty 
relation. It is therefore intuitive to expect, that out of all pure states the coherent states are the most distant from 
Psf. In the other extreme case, both zeros lay at the antipodal points, x = tt, which corresponds to /5|o) , and R^ — i?2- 
In this symmetrical case, the Husimi distribution is as delocalized as possible, and we conjecture that for every pure 
state its distance to is larger than or equal to R2. 

Thus, considering the Monge distance from p*, one gets a foliation of the space of pure states V = CP^, as shown in 
Fig. 6b. As a running parameter we may take the angle X: which describes a pure state in the stellar representation. 
This foliation is singular, since the topology of the leaves depends to the angle. The angle x = represents the 
sphere of coherent states (~ SO{i)/ SO{2)), intermediate angle represents a generic 3D manifold RP'^/Z2 (a 
desymmetrized Stiefel manifold ^ S'0(3)/Z2) of pure states of the same Xi while the limiting value x = denotes 
the M.P'^ SO{i)/0{?>)) manifold of states rotationally equivalent to p\Qy Similar foliations of V discussed in other 
context may be found in Bacry and in a recent paper by Barros e Sa p8| . For comparison in Fig. 6a we present 
the foliation of as regards the Monge distance from p^. 

a) N = 2 b) N = 3 




tates |1,0 >) 



FIG. 6. Foliation of the sphere along the Greenwich meridian (a), foliation of the four- dimensional space of the A'' = 3 pure 
states along the angle x between both zeros of the corresponding Husimi distribution. The poles correspond to the distinguished 
2D submanifolds of CP^; the manifold of coherent states and the manifold of states equivalent to |1, 0). 
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FIG. 7. Sketch of the structure of the space of the mixed states A4 for j — 1 induced by the Monge metrics. Manifold of 
the coherent states C = is represented by a circle of radius Ri = 37r/16 ~ 0.589 centred at p,. Pure states isomorphic to 
P|o) are situated R2 — 1/6 ~ 0.166 from p,. Dots at smaller circles represent the positions of zeros zi and Z2, which determine 
each pure state in the stellar representation. 

Since it is hardly possible to provide a plot of A4 revealing all details of this non-trivial, 8-dimensional space of 
mixed states, we can not expect too much from Fig. 7, which should be treated with a pinch of salt. As it was 
discussed in Sect. II, from the point of view of the standard metrics, the four-dimensional manifold of the pure states 
V is contained in the sphere 5^ centered at p*. For the Monge metric one has Ri > R2, so we suggest to illustrate 
A4 as an 8-dimensional full 'hyper-ellipsoid'. Pure states p+ and p_ occupy its poles along the longest axis. The 
dashed vertical ellipse represents the space of all coherent states C, which forms the sphere S"^ of radius Ri. Solid 
horizontal ellipse represents these pure states, which are closest to p*. This subspace may be obtained from |0)(0| by 
a three-dimensional rotation of coordinates; topologically it is a real projective space MP^. Although both ellipses do 
cross in the picture, both manifolds do not have any common points, what is easily possible in the four-dimensional 
space v. To simplify the identification of single pure states we added in the picture small circles with two dark dots, 
which indicate their stellar representations. In general, the states represented by points inside the hyper-ellipsoid are 
mixed. However, since M is only a part of the hyper- ellipsoid, not all points inside this figure do represent existing 
mixed states. 
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FIG. 8. As in Fig. 4 for j = 2. 



E. The cases: j=3/2 and j=2 

For j = 3/2 (N = 4) the results read in a simplified notation D^r^^j^ = ^^-{j2^-3/2 ~ 57r/32, D^^^f[^2 ~ 97r/64, 
£)^^2=;*^ = D^J^l = 297r/128, and £'(72^,^^ = D^_^=^l ~ 0.2737. For j = 2 {N = 5) one obtains i'^^^^^ = 

Di'[^_'l = 357r/256, 2?^;^=^^ = D^.'^'s,'^ = 157r/128, d(;!='^) = iPL^=^) = 65^/256, i^^^I^'^ = oi^i^'^ ^ 0.3909, and 

Dffl~^^ = 29/120. Fig. 5c, d presents a schematic map of these states. Although the results are analytical, we give 
their numerical approximations, which give some flavour of the geometric structure induced by the Monge metric. 

F. Monge distances for an arbitrary j 

1. Eigenstates o/J^ 

Using the formula for distribution functions one may express the distances between neighbouring eigenstates 
of Jz for an arbitrary j by the following formula 



DM{\j,m), \j,m- 1)) = (';)2-^- (4.14) 

for m = — J ■ + 1, . . . , j, where = 2j + 1 and n = j + m=l,...,N— 1 (for the proof see Appendix El). This leads 
to the following asymptotic formula 

DM{\j,m),\j,m-l})r^^=^== (4.15) 

\Jn{N — nj 
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valid for large j, where n is defined above. It is easy to show that that for each j all the eigenstates of Jz are located 
on one metric line. Hence we get 



Dm{p+,P-) = DM{\j,j), \j, -j)) = XI ^M{\j,m), \j, 

m = -j+1 



m ■ 



1)) 



2N 
N 



-2N 



(4.16) 



2. Distance from p* 



According to Property B the distance of each coherent state from p* is the same and equal to Ri 
This quantity may be found explicitly for an arbitrary A'' = 2ji + 1: 



1 TT 

Dm{p+,P*) = -Dm{p+,P-) = - 



1 - 



2N 
N 



21- 



-2N 



Dm{p+,P*)- 



(4.17) 



which is asymptotically (for large N) equal to 7r/2 — ^J-k/N . Such a quantity is shown in Fig. 8 (for j = 2) as the 
area between two corresponding distribution functions. Observe that in comparison with Fig. 4 the coherent states 
are more localised, and the area between steeper distribution functions is larger. In the classical limit TV ^ oo we 
arrive at Dm{p+, p*) 7r/2 and Dm{p+, P-) — * tt. The latter result has a simple interpretation: in this limit the 
coherent states become infinitely sharp and the Monge plan consists in the rotation of the sphere by the angle tt. The 
three points p+, , and p_ form another metric line, which for TV > 2 is different from the metric line generated by 
the eigenstates of Jz- Thus, for N > 2, the metric induced by the Monge distance is not 'flat'. Moreover, for j € N 
we have 



(4.18) 



which tends to 7r/2 — 1 in the semiclassical limit j ^ oo (for the proof see Appendix E2). It is well known that this 
convergence is very slow. 

3. Coherent states 

Now, let us consider two coherent states <p) and ip'). It follows from the rotational invariance of the Monge 

metric (Property B) that their distance depends only on S - the angle between two vectors on S'^ representing these 
coherent states, and is equal to the Monge distance between two coherent states lying on the Greenwich meridian 
p+ = |0, 0)(0,0| and ps :~ |E!,0)(E!,0| (for N = 2 the latter corresponds to the state labelled in Fig. 3 by pe). We 
denote this distance by C(S, j) := Dm{p+, Pe)- Using Propositions 6 we obtain the following formula for this quantity 
(for the proof see Appendix E3): 



C{E,j)=nsm{E/2)Wj{sm^ (S/2)) 



where Wj is a polynomial of the form 



4J+1 

0<.u,v 
u-\-v<j 



(4.19) 



(4.20) 



The symmetric coefficients Sj^u^y are given by 



3,u,v - -2{u + v)- ly.ulvl {u + v + 1)!4»+'' 
and the asymmetric coefficients by 



(4.21) 
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oo /2s\ 

Auv-= V 7 — • (4.22) 

s—v+l ^ ' 



Note that A^^v can be also written as a finite sum 

22M+1 " C^s 



A. 



a")(2u+l) g(.+ lV.)4^- ^^-'^^ 



The rank of Wj is [j — 1/2J, i.e., the largest integer less than or equal to j — 1/2. We have W1/2 (a;) = \ (and so 
C(S, 1/2) = (7r/4) sin(S/2) ~ i?ffs(p+, Ps)), VFi {x)^% W^/^ {x) = ^ (57 + x), 1^2 [x) = (25 + x), etc. 

One can show that all the coefficients of the polynomial Wj are positive. This leads to the following simple lower 
and upper bounds 

TT Cj sin (S/2) < C (S, j) < n Dj sin (S/2) , (4.24) 

where 

Q := 1^,(0) - 3^2 ([3/2, 1/2 - J, 1 - j1 , [2, 2] , 1) , (4.25) 

and 

7 f27 + 1) 

D, ■.= W,il) = ''J^^^ ' (2- 3^2([l,l/2-j,l-j-],[2,3/2],l)- 3^2([l,l/2-j,l-j],[2,2],l)) . (4.26) 

(here 3J^2 stands for the generalised hypergeometric function). Note that Cj 2/7r and Dj ^ 1 in the semiclassical 
limit J — > 00. For two infinitesimally close coherent states the angle S ~ and we get 

C(S,j)~|qS, (4.27) 

with 

C(S,j)^S (j^oo), (4.28) 

which agrees with Property A. 

4- Chaotic states 

In the stellar representation the coherent states are represented by iV — 1 zeros merging together at the antipodal 
point on the sphere. These quantum states are rather exceptional; a typical state has all zeros distributed all over 
the sphere. It is known [M that for the so-called chaotic states (eigenstates of Floquet operator corresponding to 
classically chaotic systems) the distribution of zeros is almost uniform in the phase space. Such states are entirely 
delocalized and their Husimi distribution is close, in a sense of the Li metric, to the uniform Husimi distribution Hp^ 
corresponding to the maximally mixed state p*. One can therefore expect (applying Proposition 2) that the Monge 
distance between these chaotic pure states pc and is small. We conjecture that the mean value of the Monge 
distance Dm{Pc, P*) of randomly picked chaotic state pc from p^ tends to in the semiclassical limit j 00. 

G. Correspondence to the Wehrl entropy and the Lieb conjecture 

In order to describe the phase space structure of any quantu m sta te p it is useful pQ| to define the Wehrl entropy 



Sp as the Boltzmann-Gibbs entropy of the Husimi distribution (4.6) 



Sp = - Hp{9, ifi) HHp{9, ^)]dpid, ^). (4.29) 
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It was conjectured by Lieb [|6l| that this quantity is minimal for coherent states, which are as locahsed on the sphere 
as it is aUowed by the Heisenberg uncertainty relation. For partial results in the direction to prove this conjecture 
see p2[-|65|. The minimum of entropy reads S'min = {N — 1)/N — InA^, where the logarithmic term is due to the 
normalisation of the Husimi distribution. It was also conjectured |]63|] that the states with possibly regular distribution 
of zeros on the sphere, which is easy to specify for Pythagorean numbers N = 2,4,6,8,12,20, are characterised by 
the largest possible Wehrl entropy among all pure states. 

Let us emphasise that for iV >> 1 the states exhibiting small Wehrl entropy comparable to S'min are not typical. 
In the stellar representation coherent states correspond to the coalescence of all iV — 1 zeros of Husimi distribution 
in one point. In a typical situation all zeros are distributed uniformly on the sphere, and the Wehrl entropy of such 
delocalized pure states is large. Averaging over the natural Haar measure on the space of pure states V one may 
compute the mean Wehrl entropy (S) for the A^— dimensional states. In slightly different context such integration was 
performed in ||66|-|6^ leading to 

(S)^(^) =-ln7V + v|/(7V+l)_vI;(2) , (4.30) 

where \I' denotes the digamma function, which for natural arguments k < n satisfies 'ii{n) — 5'(fc) = J2m=k I/"*- 
the classical limit N ^ oo the mean entropy tends to 7 — 1 ^ —0.42278 (7 is the Euler constant), which is close to 
the maximal possible Wehrl entropy Sp^ = 0. 

The Wehrl entropy does not induce a metric in the space of quantum states. However, it describes the localisation 
of a quantum state in the classical phase space and has some properties similar to the Monge distance of a given 
state p to the maximally mixed state [ p9[ . In view of our results on the Monge distance, we advance analogous 
conjectures, concerning the set of pure states V belonging to the A^— dimensional Hilbert space. 



Conjecture 1. In the sense of M onge metric the coherent states are pure states most distant from p*. This 



maximal distance Ri is given by ( 4.17) and tends to tt/2 for N ^ 00 



Conjecture 2. Pure states which maximise the Wehrl entropy are the most close to p* in the sense of Monge 
metric. This minimal distance R2 is equal to 1/6 for N — 3 and tends to for A^ — > 00. 



In the analogy to the properties of the Wehrl entropy and formula ( 1.3C| ), one can expect that the mean distance 



(R) = (£'m(|V')(^I) P*)) averaged over the natural measure on the manifold of pure states V, is close to the minimal 
distance R2 and, for large A, is much smaller than the maximal distance In other words, the coherent states, 
distinguished by the fact of being situated in M as far from as possible, are not generic. This observation is not 
surprising, since C ~ CP^ while V ^ CP^^^, but is not captured using any standard metrics in the space Ai of 
mixed quantum states. 



V. COMPARISON OF MONGE AND STANDARD DISTANCES 



Results obtained for distances between several pairs of mixed states are summarised in Table 2. Calculation of the 
trace, Hilbert-Schmidt and Bures distances are performed directly from the definitions provided in the Sect. I. 

Note that the geometry of the space Ai is well understood for A^ = 2. In the sense of the trace and the Hilbert- 
Schmidt metrics the set of all mixed states has then the property of a ball contained inside the Bl och sphere: the 
states p+, p* and p_ form a metric line. The same statement is true for the Monge metric (see formula 4.11). However, 
for the Bures metric the situation is different. As it was shown by Hiibner the set M has in this case the structure of 
a half of a 3- sphere ||], so p+,p* and p_ form an isosceles triangle. However, the state p*, located at the pole of S^, 
is equally distant (with respect to Bures metric) from all the pure states V, which occupy the 'hyper-equator' ~ S^. 



A. Geometry of quantum states for large N 



The data collected in Tab. 2 allow us to emphasise important differences between the geometry induced by the 
standard distances and the Monge distance. From the points of view of all three of the standard metrics, the distance 
R between p* and any pure state is constant. Therefore, in these standard geometries, the coherent states are not 
distinguished in any sense in V. 
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On the other hand, a 'semiclassical' geometry, induced by the Monge metric in the (iV^ — 1)— dimensional space Ai, 
distinguishes the space of coherent states C ~ 5*^. Their Monge distance (Ri) from the centre is maximal. If we 
try to visualise the {2N — 2)— dimensional space of pure states V ~ CP^~^ as a 'hyper-ellipsoid', the coherent states 
form the 'largest circle', represented by the dashed ellipse in Fig. 9. There exists also a multidimensional subspace of 
V, consisting of delocalized pure states pc, with zeros of the corresponding Husimi function distributed uniformly on 
the sphere. Such states are situated close to with respect to the Monge metric. In the classical limit N oo, their 
distance from (i?2) is arbitrary small, so the manifold V, almost touches the maximally mixed state p*. In this 
case, we might think of as of a full (iV^ — 1)— dimensional disk of radius Ri ~ 7r/2 centred at p*, with coherent 
states at its edge and the pure states on its surface. Since it is rather flat, and contains a lot of its 'mass' close to its 
centre, it resembles, in a sense, the Galaxy. 



States 



D 



HS 



D 



Bure 



D 



Monge 



V2 



V2 



7r(l 



nil/2- f^)2-^^)^vr/2-v/V]V 



{P+,P*) = {P-,P*) 
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line 
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line 
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1 
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A^-dim simplex 
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{P+,Pa) 



2(1 - a) 



V2{1 



2vr(l- f^)2i-2A')(l-a) 



^2(1- vT^y 



{P-,Pa) 



2a 



V2a 



2^(1- f^)2i"2A')„ 



P+,Pa,P- 



line 



line 



A 



line 



(p+,Ph) 



iV = 2 
iV > 2 



2sin(S/2) 
2^1-cos4j(S/2) 



V2sin(S/2) 
^2-2cos4J(S/2) 



^2-2cos(S/2) 
^2-2cos2j(S/2) 



(7r/4)sin(S/2) 
C(S, j) = TTsin (S/2) W, (sin^ (S/2)) 



(P-,Ph) 



N ^2 
N>2 



2cos(S/2) 



2./l-sin4^'(S/2) 



\/2cos(S/2) 
2-2sin'^J'(S/2) 



y/2 - 2sin(S/2) 



2-2sin2J(S/2) 



(7r/4) cos(S/2) 
C(S, j) = ^ cos (S/2) J, (cos2 (S/2)) 



Table 2. Standard distances (trace, Hilbert-Schmidt and Bures) versus the Monge distance for various quantum 
states in TV = 2j + 1 dimensions - pure states: the coherent states ps — |S)(S|, p^, p_, and the eigenstates |m) of 
Jz, and mixed states: pa, defined in Sect. IVB, and the maximally mixed state p*. For the Monge distance we give 
semiclassical asymptotics (N oo). The polynomials Wj{x) are given by formula (4.2C). 



B. Dynamical properties 



As mentioned in Sect. II the standard distances are preserved by the unitary dynamics (see formula ( 2.11 )). Anal- 
ogous relation is true for the Monge distance only for some special cases, e.g., for simple rotations U ~ exp(iaJz) 
which preserve the coherence. In general, however, the Monge distance is not conserved 

Dm{pi,P2)^Dm{p[,p'2)- (5.1) 

Vaguely speaking, during the rotation of the 'hyper-ellipsoid', depicted in Fig. 9, a kind of contraction occurs, so 
the Monge distance changes during the unitary time evolution (and hence is not a monotone metric). Since in the 
classical limit the distance between coherent states tends to the classical distance on the sphere, we suggested Q 
to study the time evolution of the Monge distance Dj\f (t) between two neighbouring coherent states. The quantity 
A(t) = lim£)j^j(o)^o(lii[I?A/(t)/I?A/(0)])/i characterises indeed the stability of the quantum system. To get a closer 
analogy with the classical Lyapunov exponent one should then perform the limit t — s- oo. However, for longer times, 



22 



both vector coherent states become delocahzed (under the assumption of a generic evolution operator U), and their 
distance to becomes small. Therefore, after some time tr, the distance DM{t) starts to decrease, so instead of 
analysing Hmj_+oo A(t) (which always tends to zero), one needs to relay on a finite times quantity X{tr) pl| , [70| . 




FIG. 9. Sketch of the space M in the semiclassical regime j >> 1. In the limit j ^ oo the larger radius Ri tends to 7r/2, 
while the smaller R2 — » 0. Dots and small circles show the corresponding pure states in the stellar representation. 



C. Delocalisation and decoherence 



As mentioned above, the localisation of a given pure state |0) in the classical phase space is reflected by its large 
Monge distance from p*. In an analogous way one may characterise the properties of a given Hamiltonian H or a 
unitary Floquet operator F by the mean distance of its eigenstates \vi), i — I,. . .,N, from the maximally mixed 
state. Such a quantity, 7 := J2i=i DM{\vi){vi\, p^)]/N, indicates the average localisation of the eigenstates, relevant 
to distinguish between integrable and chaotic quantum dynamics f7l|| . It might be thus interesting to find unitary 
operators -Fmin and -Fmax, for which the mean distance 7 achieves the smallest (the largest) value. 

Physical systems coupled to the environment suffer decoherence. The density matrix of a given system tends to be 
diagonal in the eigenbasis of the Hamiltonian Hj, which describes the interaction with the environment | |7^ . In the 
simplest case, N ~ 2, the decoherence may be visualised as an orthogonal projection into an axis determined by Hj. 
For example, if Hj is proportional to Jz, it is just the z axis, which joints both poles of the Bloch sphere. 

In the general case of arbitrary TV, there exists an (N — 1)— dimensional simplex X of density matrices diagonal in 
the eigenbasis of Hj. Decoherence consists thus in projecting of the initial state into T. In a generic case of a typical 
interaction the eigenstates of Hj are delocalized and their Monge distance from is small. On the other hand, the 
typical coherent states are located far away from X, in the sense of the Monge metric. One can therefore expect, that 
the Monge distance of a given quantum state from X contains the information concerning the speed of decoherence. 
It is known that among all pure states the decoherence of the coherent states is the slowest fl^ . 

Moreover, the speed of decoherence of a Schrodinger cat-like pure state, localised at two different classical points x 
and x' , in a generic case depends on their distance in the classical phase space. Consider now a coherent superposition 
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l^/i) = + \f3))/y/2 of arbitrary two quantum pure states. The Monge distance between them, DM(|a), might 
be thus used to characterise the speed of the decoherence of the cat-like state 



VI. SIMPLIFIED MONGE DISTANCE BETWEEN PURE STATES 



A. Definition 



With help of the stellar representation [ fS3| , p6y55| ] we may link any pure state \lp) of the A^-dimensional Hilbert space 
to a singular distribution fip{x) containing [N — 1) delta peaks placed in the zeros Xi of the corresponding Husimi 
function H\^^^\(x), where x £ 5*^, 



AT-l 



IV'>-^(a;):-^^^(5(a;-xO. (6.1) 

i=l 

The zeros Xi may be degenerated. For any coherent state all {N — 1) zeros cluster at the antipodal point, so \a) is 

represented by fa{x) — 6(x — a). 

The simplified Monge distanc e bet ween any pure states \(p) and is defined as the Monge distance (3.3) between 



the corresponding distributions (6.1 



D,M{\v),m -^DMiUJi,). (6.2) 

It may be also called discrete Monge distance, since it corresponds to a discrete Monge problem, which may be 
effectively evaluated numerically by means of the algorithms of linear programming EH| . This contrasts the original 



definition (3.11), for which one needs to solve the two dimensional Monge problem for continuous Husimi distributions. 

Clearly, in the space of pure quantum states both Monge distances are related. This fact becomes more transparent, 
if one realises that ( |6.2[ ) is equal to the Monge distance between the related distributions f^p := i^iLi^ ^Xi]/{N — 1) 
and := Etii^ ^yi]/{-^ ~ 1)j where j/i, . . . , un-i are zeros of the Husimi distribution H^^^^^, and points Xi and 
Xi (resp. Hi and yi) are antipodal on the sphere (note that the bar does not denote here the complex conjugation). 
The distributions and may be considered as a discrete, (N — l)-points approximation of the continuovis Husimi 
distributions H^^^^^^ and H^^^^^. 

Since any coherent state is represented by a single Dirac delta, faix) — 6{x — a), the semiclassical condition ( |l.8| ) 
is exact for any dimension N 

^sA/(hl),|?72))=d(?7l,?/2). (6.3) 



Thus for N — 2 the discrete Monge distance, Dsm, is equal to the Fubini-Study distance (L7), (in this case, the 



Riemannian dista nce d on the sphere), while the continuous Monge metric, -Dj\/, is proportional to the Hilbert- 



Schmidt distance (1.2), (in this case, the Euclidean distance along the cord inside the sphere). At small distances 



both geometries coincide (the 'flat earth' approximation). 



B. Eigenstates of Jz 

In stellar representation the state |j, to) is described by j + to zeros at the south pole and j — m zeros at the 
north pole. Thus the distribution fj^m{x) consists of two delta peaks, tt apart, and it is straightforward to obtain the 
following general result 

^.m(|j, m),|j,m')) = ^|m-m'|. (6.4) 

In particular DsM{\j^i)^ |i, = — d{a, a). The zeros of the Husimi function of the eigenstates of the operators 
Jy and Jx are located at the equator at the distance 7r/2 from both poles. Thus 

TT 

DsAii\j,m)^,\j,m')y) ^ DsM{\j,m)^,\j,m")x) ^ - (6.5) 
for any choice of quantmn numbers m,m' and to". 
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C. Random chaotic states 



Eigenstates of classically chaotic dynamical systems may be described by random pure states |7l|. Expansion 
coefficients of a chaotic state \'ipc) in an arbitrary basis may be given by a vector of a random unitary matrix, distributed 
according to the Haar measure on U{N). Zeros of the corresponding Husimi representation are distributed uniformly 
on the entire sphere |5^, (with the correlations between them given by Hannay ||7^). This fact allows one to compute 
the average distance of a random state to any coherent state 

D,M{\a)M)^]^j^ SsinSdS=|. (6.6) 

In a similar way we get the average distance to the eigenstates of Jz 

TTITT 

iVc)) = xsinx + cosx, where x = (6-7) 

It admits the smallest value equal to unity for m = x = 0, while the largest value is obtained for m ~ ±j, for which 
the above formula reduces to (|6.6| ). 

Let us divide the sphere into cells of diameter proportional to a/ZV- Consider two different uncorrelated random 
states lipc) and 10c)- Uniform distribution of zeros implies that there will be on average one zero in each cell and 
the distance between the corresponding zeros of both states is of order of ViV- Thus their simplified Monge distance 
vanish in the semiclassical limit, 

D,Mm,m)^^^-N-^/^^0 (iV^oo). (6.8) 



Thus in the space of pure quantum states the simplified, discrete Monge metric Dsm displays several features of the 
original, continuous Monge metric Dm- 



VII. CONCLUDING REMARKS 



In this paper we analysed the properties of the set of all mixed states constructed of the pure states belonging to 
the A^— dimensional Hilbert space. The structure of this set is highly non-trivial due to the existence of the density 
matrices with degenerated spectra. Each spectrum may be represented by a point in the (A^ — 1)— dimensional simplex. 
In a generic case of a nondegenerate spectrum (point located in the interior of the simplex) this set has a structure 
of [U{N)/{U{1))^] X Gn- However, there exist all together 2^^^ parts of the asymmetric simplex of eigenvalues, all 
but one corresponding to its boundaries. These boundary points, representing various kinds of the degeneration of 
the spectrum, lead to a different local structure of the set of mixed states. 

Standard metrics in the space of quantum states are not related with the metric structure of the corresponding 
classical phase space. To establish such a link we used vector coherent states, localised in a given region of the sphere, 
which plays a role of the classical phase space. Each quantum state may be then uniquely represented by its Husimi 
distribution, which carries the information concerning its localisation in the classical phase space. We proposed 
to measure the distance between two arbitrary quantum states by the Monge distance between the corresponding 
Husimi distributions. Therefore, to compute this distance, one has to solve the Monge problem on the sphere. Thus, 
unexpectedly, a motivation stemming from quantum mechanics leads us close to the original Monge problem of 
transporting soil on the Earth surface. Even if the exact solution of the Monge problem is not accessible we can use 
either lower and upper bounds for the Monge distance (definition (3.2), Propositions 2 and 4), or numerical algorithms 
based on the idea of approximation of continuous distributions by discrete ones. These techniques lead to general 
methods of computing the Monge distance on the sphere (Propositions 5 and 6), as well as to concrete results we 
obtained in this paper (Sect. IVC,D,E,F and Sect. VA). 

The Monge distance induces a non-trivial geometry in the space of mixed quantum states. For A^ = 2 it is consistent 
with the geometry of the Bloch ball induced by the Hilbert-Schmidt or the trace distance. For larger A^ it distinguishes 
the coherent states, which are as localised in the phase space, as it is allowed by the Heisenberg uncertainty principle. 
These states, laying far away from the most mixed state p*, are not typical. The vast majority of pure quantum states 
are localised in vicinity of in the sense of the Monge metric. The Hilbert-Schmidt distance between a given state 
p and may be used to measure its degree of mixing. On the other hand, the Monge distance Dm{p,P*) provides 
information concerning the localisation of the state p in the classical phase space. 
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A similar geometry in the space of pure quantum states is induced by the simphfied Monge metric d^M- It is defined 
by the Monge distance between the {N — l)-points discrete approximations to the Husimi representation generated 
by the stellar representation of pure states. This version of the Monge distance may be easily evaluated numerically 
by means of the algorithms of linear programming |p5|] . Therefore it might be used to study the divergence of initially 
closed pure states subjected to unitary dynamics and to define a quantum analogue of the classical Lyapunov exponent 
[^,0. Moreover, this metric may be useful in an attempt to prove the Lieb conjecture: it suffices to show that for 
any pure state the Wehrl entropy decreases along the line joining with the closest coherent state. 

In contrast with the standard distances, the both Monge distances are not invariant under an arbitrary unitary 
transformation. This resembles the classical situation, where two points in the phase space may drift away under 
the action of a given Hamiltonian system. In a sense, the Monge distance in the space of quantum states enjoys 
some classical properties. Several classical quantities emerge in the description of quantum systems. We believe, 
accordingly, that the concept of the Monge distance between quantum states might be useful to elucidate various 
aspects of the quantum-classical correspondence. 

VIII. ACKNOWLEDGEMENTS 

K. Z. would like to thank I. Bengtsson for fruitful discussions and a hospitality in Stockholm. It is a pleasure to 
thank H. Wiedemann for a collaboration at the early stage of this project and to S. Cynk, P. Garbaczewski, Z. Pogoda, 
and P. Slater for helpful comments. Travel grant by the European Science Foundation under the programme Quantum 
Information (K.Z) and financial support by Polish KBN grant no 2 P03B 07219 are gratefully acknowledged. 



APPENDIX A: PROOF OF PROPOSITION 2 

Applying Proposition 1 we see that it suffices to prove the inequality 



/ 

Jn 



f{x){Qi{x) - Q2{x))dm{x) 



< {A/2)DlAQi,Q2) 



(Al) 



for every weak contraction / : ^ M. For such / we see at once that (max / — min /) < A. Let us consider a function 
g : H. ^ M. defined by the formula g{x) = f{x) — min / — A/2 for x E ft. Clearly j^l < A/2. Finally, we get 



f{x){Qi{x) - Q2{x))dm{x) 



'j{x){Qi{x) - Q2{x))dm{x) 



< {A/2)Dl,{Qi,Q2) 



(A2) 



which completes the proof. 



APPENDIX B: PROOF OF PROPOSITION 4 



Let pi, p2 G -M- We denote the set of all contractions (c— Lipschitzian functions with c < 1) f : fl 
We have 

Dm{pi,P2) = Dm {Hp^,Hp,J = 



by Lipi 



max 



max 

feLipi 



max 

feLipi 



tr 



f{v){v\{Pi - P2)\v)dm{T]) 
f{v)\v){v\dm{v) {Pi - P2) 



max •|| ti A{pi — p2)\ '■ A — Hermitian and A = J f {Ti)\ri) {ri\dm{ri) for some / G iipij- = 



max I ii A(pi — 02)! 

L{A)<1 ' 



(Bl) 
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APPENDIX C: PROOF OF PROPOSITION 5 



We start from two simple lemmas on weak contractions on the sphere. In the sequel, d denotes the Riemannian 
metric on S^. 

Lemma 1. Let / : [0, tt] — + R be a weak contraction. Define / : 5^ — * R by the formula 

mip) = fi^) (ci) 

for {■&, (p) E S"^. Then / is a weak contraction. 

Proof of Lemma 1. Let {'di,(pi),{'d2,(p2) G S^. Applying spherical triangle inequality we get 

IM, (^l) - /(l?2, ^2)1 = Ifi^l) - fi^2)\ <Wl-^2\< (Pl), (^2,^2)). 

Lemma 2. Let G : 5^ ^ R be a weak contraction. Define G : [0, tt] ^ M by the formula 



2tt 



for £ [0,7r]. Then G is a weak contraction. 

Proof of Lemma 2. Let i?i,i?2 e [0,7r]. We have (^), (i?2, <^)) = - ^2]- Hence |G(i9i) - 0(^2) 

{l/2n)\ Jo "(G(^i, ^) - G{d2,V^))dv\ < (1/2^) ^ dm^v)^ {^2,^))d^ < \A - ^2\. 



(C2) 



max 
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JO 
7r / p27v 



Proof of formula (4.7). 

It follows from Proposition 1 that 

Dm{pi,P2) = Dm {Hp^,Hp^] 

f{{d, ^)){Hp, (^) - Hp, (dj) sin d/4ndd^ d^ 
From the Salvemini formula (3.4), Proposition 1, the above Lemma 1, and formula ( |C3| ) we deduce that 

\FpM^FpM\di} = 
fi^){{hpM~hpM)di9 = 
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V5))(i?p, (i?) ~ iJp, (^?)) sin d dd] d^ 
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Dm{pi,P2) ■ 



(C4) 



On the other hand, applying formula (C3), the above Lemma 2, Proposition 1, and the Salvemini formula (|3.4[) we 



get 



max 

GeLipiiS^} 



Dm{pi,P2) 

G{{d,ip)){Hp, (i9) - i^p,(^9))sin^?/47^^^19 ) d(p 



max 

GeLipi(S2) 
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geLjpi([0,7r]) 
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G(i9) ((-ffi(i9) - i?2(i?))sini?di5 
5(7?) (iJi(i9) sini9- iJ2(i?)sini?) di? 

\Fi{d)-F2{d)\d^, 



(C5) 



which establishes the formula. 
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APPENDIX D: PROOF OF PROPOSITION 6 



Put 



Dm{pi,P2) ■■= 77- 



27r /.TT 



Jo 



{Flit,!}) ^ F2(t,i))) dt dd. 



Let (/? € [0, 27r] and i = 1,2. Integrating by parts we get 



Hence 



(Dl) 



(D2) 



(Fi (t, (^) - F2 [t, ^)) dt = / i-t) - {Hp, [t, if) - Hp, {t, if)) sint dt + t (Fi {t, ip) - F2 {t, ip))\'^l . (D3) 
Jo ^ 



According to assumption (1) the last term is equal to 0. Thus, applying Proposition 1 to / : S*^ 
/ {t, If) = -t, for t e [0,7:], If e [0, 2n], we get 



given by 



27r /"TT 



JO 



Dm{pi,P2)=I I {-t) -{Hp, {t,(p) - Hp, {t,ip))smt dt dip < Dm {pi,P2) 



(D4) 



On the other hand consider the following transformation of the density Hp, into Hp, : we transport the 'mass' along 
each meridian separately (it is feasible due to assumption (1)) and then we join all the transformations together. 
Applying the Salvemini formula (3.4) to each meridian {tp G [0, 27r]), averaging the results over [0, 27r], and finally 
using assumption (2) we get 



Dm (pi,P2) < ^ 



27r /•TT 



Jo 



\Fi {t, (p) - F2 {t, if) I dtd^^ Dm (Pi, P2) 



(D5) 



which completes the proof. 



APPENDIX E: DERIVATION OF THE MONGE DISTANCES FOR SOME INTERESTING CASES 



1. Derivation of formula (4.14) 



Let j e N and m = + We put = 2j + 1 and n — j + m. Applying Proposition 5 and then the 

substitution u = cos^ (i^/2) we get 



DmHJ, m), \j, m - 1)) = /J" G (n, N, cos^ (i9/2)) - G (n - 1, iV, cos^ {d/2)) (sini9) /2dd\ dy 



= lo iLs- 1 G {n, N,u)-Gin~ 1, N, u) du 
where G {n,N,u) := 7V(^,;^)u" ( 1 - u)'^"^""^ for u e [0, 1]. Using the identity 



dy, 



(El) 



/ 



(G (n - 1, iV, w) - G (n, N, u)) du 



N' 



■{l^u) 



(N-n) 



(E2) 



we obtain 
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DMi\j,m), |j,m- 1)) 



N\T{n+l/2)r{N -n + 1/2) 
n j fjN) 

2 (iV - n)\ /2n\ _2jv 



N — n J yn 



as desired. 



(E3) 



2. Derivation of formula (4.18) 

For j e N we put Dj :— Dm(| j, 0), p*). From Proposition 5 and formula ( [4.7D we deduce that 



^li.o>O-,o|('0)(^'0 - / hp^{il))dil) 
Jo 



(G (cos^ (VV2)) - 1) - sinV'dV 



M 



(E4) 



where G (u) {2j + 1) (^j')w^ ( 1 — m)-' for m e [0, 1]. Applying the substitutions u — cos^ (V'/2) and y — -0/2 and 
using the symmetry arguments yields 



D,^2 







7r/2 

(G(li) - 

cos'^(iJ/2) 

= 2 / / G(u) dudy-l 

•/ siii^ y 

7r/4 /"cos^ y 



sm y 



( 1 — u) dudy — 1 . 



(E5) 



Set cj {u) := 2 (2j + 1) (2/) / u^' ( 1 - u)^' du for u e [0, 1], j G N. Then Dj = J^' {cj (cos^ y) - (sin^ y)) dy - 1. 
Integrating by parts we get Cj (u) ~ 2( j') (2u — 1) (1 — u)-' + Cj_i (m), and so £>j ~ (^j')2^^-' ^JpY + ^j-i- Moreover 
we can put Do = 0. Thus Dj — 2FFT^^^'^( fc^) ~ 2FFT ^^(2fcM!'' ' claimed. Applying Taylor's formula 



arcsma; = z^^^o 



— 1 (2fc--l)!! _^2fc+l 



2fe+l (2fc)!! 



we obtain I?j — )■ 7r/2 — 1 (j ^ 00). 



3. Derivation of formula (4.19) 



Let G(S,j) = Dm{p+, P'E.) for S e [0, tt] and j = 1/2, 1, ... . It follows from the rotational invariance of the Monge 
metric (Property B) that G(S,j) = D m{p i, P2), where pi = P(t_h)/2 and p2 = P{n+S)/2- To apply Proposition 6 
observe first that according to formula (4.5) we have 



Fp,(^,^) = (2j- + 1) 



1 + sin?? cos {^/2) cos If + cos?? sin (^/2) 



2j 



(E6) 



1 + sin ?? cos {^/2) cos ip — cos ?? sin (^/2) 



2i 



(E7) 
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and so Hp^{'d,(p) = Hp^{Tr — ^,(p) for <p) € S"^. Thus, applying the substitution n — d i}, we get Fi{tt,lp) = 
\ Hp^{d, sint^ d'd'^\ £ Hp,{d, ^) siTid dd = F^iiT, and Fi{t, <p) - F^it, ip) ^ Fi{tt ^ t, ^) ^ F2{tt ~ t, cp) for 
t G [0, tt] and </? G [0, 27r], which implies the assumption (1). Moreover, Hp^ tp) — Hp^ {■&, (p) > for i) £ [0, n/2] and 
ip £ [0, 27r]. From this fact and from the symmetry of the functions Fi{-^ip) — F2{-,(p) {ip £ [0, 27r]) we deduce the 
assumption (2). Hence the assumptions of Proposition 6 are fulfilled and we conclude that 

C{E, j) = J^^ £ (^J^ (^{w + zf' ~{w-zf') sin d dd^ dt d^ , (E8) 

where w 1 + siniJcos (S/2) cos(p and z := cosi? sin (S/2). Applying the identity 

, ^ ,2, , for 2, -even 

[w + z) - [w - z) ■' = i (E9) 

l2^ECo^' (2^^)w2fc^2j-2fe-l fo^ 2j - odd 

and performing the integration we get after tedious (but elementary) calculation the desired result. 
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